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The charm contribution to the atmospheric fluxes of muons and muon neutrinos may be enhanced by as much 
as a factor of 10 when one includes the contributions of D — > n,K -^leptons and folds in uncertainties in the 
charm cross section and energy distribution. In the energy range considered here, from 100 GeV to 10 TeV, the 
charm contribution is small compared to the conventional flux of muons and muon neutrinos. 



The fluxes of leptons from the decays of pions 
and muons produced by cosmic ray interactions 
in the atmosphere are known to within approx- 
imately ±20% [0 at energies ~ 1 GeV. At cos- 
mic ray energies greater than a few GeV, charm- 
anticharm pairs can be produced. The semilep- 
tonic decays of charmed mesons and baryons 
which emerge from the cosmic ray interactions 
with air are additional contributions to the atmo- 
spheric lepton fluxes. Here, we present the charm 
contribution to the atmospheric lepton fluxes in 
the energy range of 100 GeV to 10 TeV. We 
evaluate the theoretical uncertainties associated 
with charm production and decay. We show that 
the lepton flux from charm decays has a factor 
of approximately ten uncertainty, but the uncer- 
tainty has little implication for the measured at- 
mospheric muon and muon neutrino fluxes in this 
energy range. 

Atmospheric leptons from pion and kaon de- 
cays are called "conventional" leptons. Leptons 
from charm decays contribute to the "prompt" 
flux. Recently, a new calculation of the prompt 
lepton fluxes by Thunman, Ingelman and Gon- 
dolo (TIG)|| has appeared in the literature. Us- 
ing a PYTHIAU based Monte Carlo, they have 
evaluated the contributions of semileptonic de- 
cays of charmed particles to the leptonic fluxes. 
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One aspect of charm decays not included in 
other calculations is the decay chain of charm—* 
ir,K ^leptons, called here "secondary" contribu- 
tions. Each charm decay has an associated pion 
and/or kaon in the final state, which itself decays 
leptonically. 

The calculation of prompt, conventional and 
secondary lepton fluxes can be done approxi- 
mately using the Z-moment form of the cascade 
equations. Details of approximate solutions can 
be found in, for example, Ref. The general 
form for the flux of particles of type j, as a func- 
tion of energy and slant depth in the atmosphere, 
X, is 
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The quantities Xj and are the interaction 

and decay lengths, (Nj) is the average particle 
multiplicity of type j, and dn k ^j /dE, describes 
the energy distribution of particle j given its pro- 
duction by particle k with energy E k . Z- moments 
are defined by 



an -\ — <t>k{E,X) 
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The Z-moments Z k j depend on particle type and 
energy. 
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The starting point for the solution to the cas- 
cade equations is the cosmic ray flux, which we 
take as all protons with an energy behavior cj) p ~ 
E~ 2 - 7 — E~ s , as in TIG. Given the cosmic ray 
flux, and the interaction and decay moments of 
TIG in Ref. §, we can evaluate the conven- 
tional and prompt fluxes. Since our interest in 
is the "low energy" regime, 100 GeV< E < 10 
TeV, we include only c = D° 7 D° , Z? ± in our 
charm contribution. TIG have shown that these 
are dominant in this energy range. The new fea- 
ture here is to include secondary decays: we have 
additionally, for example, D — > it — > fi, so we 
need the decay moment Zjj^. As a first approx- 
imation, we rescale the decay moments for D 's 
into neutrinos to account for hadronic branching 
fractions and multiplicities. The details of the D 
decay- moment inputs can be found in Ref. || . 

Our results for the fluxes of muons (particles 
plus antiparticles), including secondary decays, 
are shown in Fig. 1 for the vertical direction. The 
solid line indicated by P (prompt) include only 
L>°, D° and D ± . The dashed lines are the TIG 
parameterization of their Monte Carlo results ||] , 
which have significant A c and Df- contributions 
at high energies. At low energies, the D's domi- 
nate. The muon neutrino fluxes are similar. The 
prompt flux of muon neutrinos equal the prompt 
flux of muons. The conventional muon neutrino 
flux is approximately a factor of 10 suppressed 
relative to the muon flux. The secondary muon 
neutrino flux is about a factor of five suppressed 
relative to the secondary muon flux. 

As can be seen from the figure, the secondary 
flux is approximately three orders of magnitude 
below the conventional one. As one goes to an- 
gles off the vertical, both the secondary and con- 
ventional fluxes increase, but in a constant ratio, 
while the prompt flux remains constant as a func- 
tion of angle. 

Uncertainties in the calculation of the charm 
cross section and energy distribution affect pre- 
dictions of both the prompt and secondary fluxes. 
We estimate that the cross section data || can 
accommodate an additional factor of two in the 
theoretical prediction of TIG [Q] . The energy dis- 
tribution in the charm production cross section 
is another uncertainty. The energy distribution 
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Figure 1. Conventional (C), prompt (P) and sec- 
ondary (S) atmospheric muon plus antimuon flux 
in the vertical direction. The dashed lines are the 
approximate formulae of Ref. [2] . 



is usually written in terms of the charm particle 
energy divided by the incident nucleon (beam) 
energy: x = E c /Eb, and in the scaling approx- 
imation, da/dx ~ (n + 1)(1 — x) n . Next-to- 
leading order (NLO) perturbation theory, when 
fit to the (1 — x) n distribution, has n ~ 6 — 9.5 
for E b = 100 - 1000 GeV, while experimental 
measurements yield n ~ 4.9 — 8.6 in the same en- 
ergy range Q. By using experimental rather than 
theoretical values for n, the charm flux can be 
enhanced by a factor of 1.5. 

Finally, the charmed meson decay moment 
used for Fig. 1 is based on a parton V-A for- 
mula. If we use phase space instead of the V-A 
formula, the Z decay moment is enhanced by a 
factor of 2.4. Taken together, these enhancements 
can increase the secondary flux by a factor of 7, 
and the prompt flux by a factor of 3. 

To estimate the effect of the high energy cross 
section on the low energy flux, we use a cross 
section which becomes 0.la pp at high energies as 
suggested by Zas et al.Q], where a pp is the total 
pp cross section d|: 



a LE + 0.1a pp y ' 

where a LE , for E = 100 - 1000 GeV, is the 
next-to-leading order (NLO) charm pair produc- 
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Figure 2. The charm- anticharm cross section ac- 
cording to Eq. (4) (dashed line). The solid line is 
the NLO charm cross section using CTEQ3 par- 
ton distribution functions with /i = m c = 1.3 
GeV and the dot-dashed line is <r vv . 



tion cross section, evaluated at factorization and 
renormalization scales fi equal to m c = 1.3 
GeV, using the CTEQ3 parton distribution func- 
tions!]]. We use a power law extrapolation for 
E > 1 TeV. The cross section of Eq. (4) is rep- 
resented by the dashed line in Fig. 2. As an ex- 
treme, we take n = 4 in da/dx ~ (1 — x) n , barely 
consistent with the low energy data. Our results 
for the muon plus antimuon fluxes are shown in 
Fig. 3. 

In conclusion, there is a factor of about ten 
uncertainty in the prediction for the flux of muons 
from the decay of charm in the energy range 100 
GeV-10 TeV. Relative to the prompt muon flux, 
the secondary decay contribution is significant, 
however, relative to the conventional flux, it is 
not. These conclusions also apply to the muon 
neutrino fluxes from charm decay. Because the 
conventional electron neutrino flux is small, the 
prompt flux is more important at 10 TeV than 
in the muon neutrino and muon case. This is a 
topic under further investigation by the present 
authors. 
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Figure 3. A comparison of the prompt (P) and 
secondary (S) muon fluxes using TIG parame- 
ters (solid line) and the fluxes computed using 
the cross section of Eq. (4) and a scaling energy 
behavior with n = 4. 
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